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Abstract 

In order to understand the nature of the lenses that generate microlensing 
events, one would like to measure their mass, distance, and velocity. Unfortunately, 
current microlensing experiments measure only one parameter of the events, the 
characteristic timescale, which is a combination of the underlying physical param- 
eters. Other methods are required to extract additional information. Parallax 
measurements using a satellite in an Earth-like orbit yield the projected velocity of 
the lens: v = v/(l — 2;), where v is the transverse velocity (speed and direction) of 
the lens relative to the Earth-source line of sight, and z is the ratio of the distances 
to the lens and the source. A measurement of v could distinguish between lenses 
belonging to the bulge and disk populations. We show that for photometric preci- 
sions of 1% to 2%, it is possible to measure the projected speed, to an accuracy 
of < 10% for over 70% of disk lenses and over 60% of bulge lenses. For measuring 
the projected velocity v, the percentages are 40% and 30%, respectively. We find 
lines of sight > 2° away from the ecliptic are preferable, and an Earth-satellite sep- 
aration in the range 0.7AU — 1.9AU is optimal. The requirements of the satellite 
for measuring the projected velocities of events towards the bulge are similar to 
those for measurements toward the LMC 
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1. Introduction 



Four ongoing microlensing searches have detected more than 100 candidate 
events, the great majority toward the Galactic bulge. The MACHO collaboration 
(Alcock et al. 1995b, 1996b) have reported 8 candidate microlensing events of 
stars in the Large Magellanic Cloud (LMC), more than expected by the known 
luminous stars of the Milky Way (Bahcall et al. 1994; Gould, Bahcall & Flynn 
1996) and LMC (Gould 1995b), but less than expected from standard spherical 
halo models (Paczynski 1986; Griest 1991). The EROS collaboration (Aubourg et 
al. 1995, Ansari et al. 1996) have reported 2 candidate events towards the LMC. 
The MACHO (Alcock et al. 1995a, 1996a), OGLE (Udalski et al. 1994), and DUO 
(Alard 1996) collaborations have reported more than 100 candidate events of stars 
towards the galactic bulge. This is also more than expected by known stars. These 
discrepancies suggest that a revision of the current models of the structure of the 
disk, bulge, and halo may be in order. 

The light curve generated by a microlensing event is described by the magni- 
fication A{x) = (x^ + 2)/x(x^ + 4)^/^, where x{t) is the position of the lens in the 
Einstein ring as a function of time, and is given by, x{t) = [uj'^it — to)^ + P'^]^^'^- 
Here to is the time of maximum magnification, (3 is the impact parameter in units 
of the Einstein ring radius, and cu = is the inverse Einstein ring crossing time. 
Thus the magnification curve is fit by three parameters, 00,(3, to. Of these three 
parameters, only cu yields any information about the lens itself; (3 and to reflect 
only the geometry of the event. The timescale is related to the Einstein ring ra- 
dius and velocity of the lens by — fe/'^- Here v — |v| is the magnitude of the 
transverse velocity of the lens along the Earth-source hne of sight. 
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where Vg, Vi, and Vq are the transverse velocities of the source, lens and observer, 
respectively, dd and dos are the distances to the lens and source, and d\s is the 
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distance between the lens and the source. The Einstein ring radius is defined, 

Te = ^^^rfos^(l - 2;), (1-2) 

C Clos 

Thus the timescale is a combination of the lens's physical characteristics: M, doi, 
and V. Because the bulge and disk of the Galaxy are characterized by different 
speed and distance distributions, separate measurements of v and dd would allow 
one to determine the component to which the lens belongs. Measurement of M 
would yield the mass spectrum of lenses in these components. Unfortunately, a 
measurement of alone provides no direct information about these characteristics. 
Some information can be obtained by considering probability distributions, but this 
method is inherently statistical in nature. Furthermore, lenses belonging to the 
Galactic bulge will produce events with time scales very similar to those produced 
by lenses belonging to the Galactic disk. As a result, measurements of te alone 
differentiate poorly between these two components. 

In order to discriminate between the bulge and disk components, one must 
consider methods of extracting more information from each event. Several meth- 
ods have been proposed, including using parallax to measure the projected velocity 
of the lens, v = v/(l — z). There are two basic methods of acquiring information 
from parallax: ground-based measurements and space-based measurements. The 
projected velocity of one lens has already been measured using ground-based par- 
allax (Alcock et al. 1995c). However, ground-based parallax is limited to those 
events for which the crossing time is large, t^k, 2 months; it is only for these events 
that the motion of the Earth induces measurable asymmetries in the light curve. 
Unfortunately, most of the events that have been detected so far have time scales 
significantly smaller than those required for ground-based parallax. 

Ground-based parallax illustrates a generic problem inherent in most of the 
methods that have been proposed for measuring the physical characteristics of 
lenses. These methods are only effective for a small fraction of events. Further- 
more, those events for which these methods are viable typically represent a very bi- 
ased subsample. Measurements of proper motions are biased towards bulge events. 
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color-shift measurements are biased towards bright lenses, and ground-based par- 
allax measurements are biased toward events with long durations. In contrast, 
satellite-based parallax measurements are effective for a majority of events, and 
are not biased towards a particular class of events. 

Two basic methods have been suggested of measuring the proper motion of the 
lens, fi = v/cIqi. These are photometric (Gould 1994; Nemiroff & Wickramasinghe 
1994; Witt & Mao 1995; Witt 1995; Loeb & Sasselov 1995; Gould & Welch 1996) 
and spectroscopic (Maoz & Gould 1994). Ideally one would like to combine mea- 
surements of both the proper motion and the parallax to determine both v and fi. 
These two parameters, along with the timescale, completely determine M, doi, and 
V (Gould 1996). (Combined measurements of the proper motion and parallax can 
be used to determine the true velocity of the lens only if the projected velocity, v, 
is determined, and not just the projected speed v.) As we discuss below, parallax 
measurements can be used to determine v in a significant fraction of events. Thus 
the mass spectrum, physical distribution, and velocity distribution of lenses could 
be determined if both the parallax and proper motion of events could be measured. 
Unfortunately, while optical interferometers may be able to measure proper mo- 
tions for many events in the relatively near future, events for which proper motion 
can currently be measured are rare. 

It is not strictly necessary, however, to measure both n and v in order to dis- 
tinguish between component populations of lenses. Han & Gould (1995) demon- 
strated that, for events towards the galactic bulge, the typical values of v are 
reasonably well separated. Thus a measurement of v alone, gathered from parallax 
information, can be used to determine the component to which the lens belongs. 
Furthermore, parallax measurements can often be used to determine not only v, 
but also the projected velocity of the lens, v. Information about v is even more 
useful in distinguishing between bulge and disk populations (Han & Gould 1995). 

Space-based parallax is the most promising method of extracting information 
about lenses. Gould (1995a) demonstrated the basic method of using information 
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obtained by a satellite in an Earth-like orbit to measure v, and also gave a rough 
estimate of the photometric precision needed. Due to parallax, the light curve of 
a particular microlensing event will be different as seen from the Earth and the 
satellite. This difference is given by the vector displacement of the lens in the 
Einstein ring, Ax = {coAt, Ap). Here At is the difference between the time of 
maximum magnification for the Earth and satellite. At — tQ — to, and A/3 is the 
difference between the impact parameters, A/3 — p' — p. The vector displacement 
is related to the Earth-sateUite separation by 

Ax=^, (1.3) 
re 

where r is the projection of the Earth-satellite separation vector R, onto the plane 
of the sky (perpendicular to the source vector s), and fe is the projected Einstein 
ring radius. 



Thus by measuring to and /3 for the Earth and the sateUite, one can determine the 
projected Einstein ring radius, and this, combined with the timescale t^, can be 
used to calculate the projected velocity of the lens, v. 

Complications arise from the fact that there is a four-fold degeneracy inherent 
in the determination of v, including a two-fold degeneracy in the projected speed, 
V. These degenerate cases arise from the fact the one does not know a priori which 
side of the lens-source hne the Earth and the satellite are on. Thus A/3 can have 
two magnitudes: A/3± = |/3' ± /3|, and two distinct signs: ±A/3± = ± |/3' ± /3|, for 
a total of four degenerate cases. This degeneracy can, in principle, be broken by 
measuring the small difference in inverse time scales, Acu — cu' — cu, caused by the 
relative velocity of the Earth and the sateUite (Gould 1995a). 

Launching a satellite for observing microlensing would be a major undertaking. 
It is often suggested that the most economical approach would be to attach the 
telescope to an already planned mission. It turns out, however, that the stability 
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requirements for observations of microlensing events are such that it would be more 
economical to launch a dedicated satellite. The primary concern, therefore, is that 
such a dedicated satellite would be sufficiently versatile and beneficial to justify the 
expense. We have already discussed what could be gained from a measurement of v 
of lenses, but it remains to be shown that a satellite with reasonable specifications 
will actually enable one to measure v, and whether these measurements could be 
made for a large class of microlensing events. There is the further consideration 
that there are two distinct fines of sight of microlensing observations, the LMC and 
the galactic bulge. The times of year when these two fines of sight can be observed 
from the ground are complementary, and one would like one satellite to be able 
to make measurements for both cases. Since the two cases are markedly different 
in the nature of events that are expected, the brightness of the source stars being 
observed, and in the geometry of the parallax measurements, it is not obvious that 
a satellite with specifications suited to LMC parallax measurements will also be 
able to make parallax measurements towards the galactic bulge. 

It is interesting to ask what the required lifetime of a parallax satellite would 
be to achieve a reasonable statistical sample of microlensing events. Gould (1995a) 
estimated that ~ 100 giant star events will be seen towards the bulge during each 
~ 6 month bulge season. To observe each event 80 times requires observing ~ 44 
events per day, or one event every ~ 30 minutes. The source stars will be typically 
/ ^ 17, thus to achieve 2% photometry requires observations of ~ 10 minutes 
on a 25 cm telescope. It is therefore possible to make measurements with 2% 
photometry for all ~ 100 events. As we discuss below, it is possible to measure the 
projected speed for ~ 70% of these events. Gould (1996) estimates that ~ 15% of 
giant events could yield measurements of the proper motion, //. Thus for a satellite 
lifetime of 3 years, we could expect to use combined parallax and proper motion 
information to measure the mass and distance of ~ 45 lenses. 

In our analysis, we choose to assume that the satellite in on an Earth-like orbit, 
separated from the Earth by a distance dsat- As we show in the Appendix, this 
orbit well approximates a realistic orbit in which the satellite has a small velocity 
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relative to the Earth at the time of launch. Such an orbit is attractive on practical 
grounds since it requires the least amount of energy to attain. Other orbits, such 
as those inclined to the ecliptic plane, or with highly eccentric orbits, require a 
much larger velocity relative to the Earth at the time of launch, and are therefore 
more expensive to achieve. An Earth-like orbit is therefore optimal for the analysis 
because it is simplest to study and is the most feasible. As we demonstrate below, 
such an orbit is also advantageous because it allows for parallax measurements 
towards the LMC and bulge equally well. 

In a earlier paper, Boutreux & Gould (1996) used Monte Carlo simulations of 
observations of microlensing events in the direction of the LMC to determine the 
conditions needed to measure the projected speed, v and the projected velocity, 
V. They found that with photometric precisions of 3% for the Earth and 4% for 
the sateUite, and an Earth-satellite separation of > 0.5AU that it was possible to 
measure v for 70% and v for 50% of LMC events. They found that larger dgat were 
preferable, and that parallax measurements were possible for a broad range of lens 
masses. 

In this paper, we extend the Monte Carlo analysis of Boutreux & Gould (1996) 
to observations of microlensing events in the direction of the bulge in order to deter- 
mine if similar satellite requirements are needed for bulge and LMC observations. 
We have determined the conditions needed to break the degeneracy for a majority 
of events, and thus to measure the projected speed, and the projected velocity, 
V, to an accuracy of < 10%. We consider the effects of the sateUite requirements, 
the photometric precision and Earth-satellite separation, as well as the mass of 
the lens and the position of the source star, in order to ensure that the parallax 
measurements will be sensitive to a broad class of events. We also consider the 
time of year the observations are made, which, as we show, is a critical parameter 
for bulge observations. As we discuss in § 2.4, it is predominantly for this reason 
that the analysis of bulge observations is substantially more complicated than that 
of LMC observations. The reader may wish to review Boutreux & Gould (1996) 
in order to become familiar with a similar problem in a simpler setting. 
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2. Measuring v 



The ability to measure v is dependent on a number of different parameters, 
several of which are tied together in intricate ways. In this section, we define 
these parameters, and use a qualitative comparison between observations towards 

LMC and observations towards the bulge to explain how these parameters affect the 
measurement of the projected velocity. In § 3 we justify these ideas mathematically. 
We quantify the ability to measure v by introducing the Degeneracy Breaking 
Fraction (DBF). The DBF can be thought of as the percentage of events for which 
the projected velocity can be measured to an accuracy of < 10% for a specified set 
of the independent parameters. 

In order to measure the projected velocity, one must be able to break the 
degeneracy between the four possible solutions. We formulate a constraint (derived 
from the difference in inverse timescales) , and discard those solutions which do not 
obey this constraint. We require that those solutions allowed by the constraint 
have fractional errors which are small, and are close enough that their difference is 
unimportant. If these requirements are met, we consider the degeneracy broken. 
We discuss the DBF quantitatively in § 2.3. 

2.1. Constraint Condition 

We distinguish the four discrete solutions by means of a constraint on the 
observables. At, Af3, and Aa;. Here Aa; = a;' — a; is the difference in inverse 
timescales between the satellite and Earth. Those solutions that are unallowed by 
the constraint at the 3a level are distinguishable. As we show below, the ability 
to break the degeneracy is complicated by the fact that the constraint coefficients 
are strong functions of the time of year. 

The constraint is (Gould 1995a) 

u\\ujAt + u_iAl3 -rAcu ^0, (2.1) 
where u is the relative Earth-satellite velocity projected onto the plane of the sky. 
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■u_i_ and arc the components of u perpendicular and parallel to r, and r = |r|. In 
order to derive expressions for -u^, u^^, and r, we define a as the angle between the 
source vector s and the south ecliptic pole, and ip as the phase of the orbit such 
that when ip — 0, the projected displacement vector r is most closely aligned with 
s. We define = 27ryr~-'^, and R — |R| as the magnitude of the Earth-satellite 
separation vector. The geometry is shown schematically in Figure 1. Using these 
definitions, we find, 

r = -sin^acos^?/;]^/^ 
uj_ = -il^R— cos a, (2.2) 



r 

u ^ ^ 



i? . 2 sin(2V>) 
— iZffiit— sm a — — - 
r 2 



Note that in Gould (1995a) the sign of W|| was incorrectly given as positive. We 
now recast the constraint equation in the form 



5]«iai = 0, (2.3) 

i=l 

where ctj are the coefficients of the constraint, and Oj are the event observables 

ai = {At,Ai3,Auj). (2.4) 
Using this, we find the simpliest form for the constraint coefficients is, 

ccj = sin^ a , Q^coso;, (1 — sin^ acos^ '0)]. (2.5) 



In order to break the degeneracy, the constraint equation (2.3) must distinguish 
between the four degenerate solutions at the 3a level. But, as we have shown 
above, the constraint coefficients are dependent on a and ip, the line of sight 
of observations and the phase of the orbit. Furthermore, from equations (1.3) and 
(2.2), we see that the vector displacement Ax is a function of a and ^. Thus, 
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ecliptic plane 




Earth 



Figure 1. The geometry of the orbits of the Earth and the sateUite. Note in particular the angle defined 
as the angle such that when = 0, the projected displacement vector r is most closely aligned with the 
line of sight. 
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because Ax = {ujAt, Af3), the event observables ai are also functions of a and 
ip. All of these contributions are obviously difficult to disentangle, and thus it is 
not a straightforward matter to describe the conditions under which the constraint 
equation (2.3) can break the degeneracy. 

2.2. Fractional Error 

We define the fractional difference as the difference between the allowed solu- 
tion and the true solution, divided by the true solution. In order that the degen- 
eracy be broken, we require that the fractional errors in the allowed solutions be 
< 10%. The behavior of the fractional errors is easier to understand than the be- 
havior of the constraint. Roughly speaking, as Ax decreases, the fractional errors 
rise. Thus, from equation (1.3) we can see that as r decreases, or as fe increases, 
the fractional errors rise. From equation (2.2), we see that the projected Earth- 
satellite separation, r, is a function of direction of the source, a, and the phase of 
the orbit, ■0- Therefore these parameters affect the fractional error, and thereby 
affect the DBF. 

2.3. DBF AND Errors 

We define the DBF as the fraction of events for which the degeneracy can be 
broken for a specified set of parameters. We consider that the degeneracy is broken 
if the difference between all allowed solutions and the true solution is less than 20% 
of the true solution, and the fractional error in Ax of each allowed solution is less 
than 10%. In principle, it could be difficult to characterize the errors in the DBF 
precisely. This is because there are two distinct types of errors that must be 
considered. The first is an intrinsic error that arises from the uncertainty in the 
values of each distinct solution. This error is approximately Gaussian distributed. 
The second error is characterized by two discrete allowed solutions, and hence is 
highly non-Gaussian. Taken together, these errors can not be described simply. 

We would like to be able to characterize our errors as 10% at the la level. 
We then consider all possible events with 10% fractional error but with multiple 
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allowed solutions. If two allowed solutions are separated by 20%, then the la error 
contour is strongly perturbed by the presence of a second solution, while the 3a 
contour is barely affected. On the other hand, if the discrete solutions are 5% 
apart, then even the la contours are essentially unaffected. Thus, at first sight, 
the threshold for the allowed separations of distinct solutions appears to depend 
sensitively on the error contour of interest. In practice, however, the DBF is very 
nearly the same for 5% and 20% thresholds. See Figure 2, below. In this paper we 
adopt 20%, which corresponds roughly to focusing on the error contours at the 2a 
level. 

2.4. LMC VS. Bulge 

Observations towards the Galactic bulge are markedly different from those 
towards the LMC. For observations towards the bulge, the results are complicated 
by the fact that the ability to break the degeneracy is a strong function of the time 
of year; this is not the case for observations towards the LMC. 

For observations towards the LMC, a ~ 0, and from equation (2.2), we find 
r — R, M|| ~ —Rfl(^, and u± — 0. These expressions are obviously independent of 
ip. Using equation (2.5), the constraint coefficients become, ctj = (O,^®, 1). Thus 
the ability to distinguish between the degenerate solutions is independent of the 
phase of the orbit, and depends only on the photometric precision, the number of 
observations per crossing time, and the Earth-satellite separation. Given a certain 
set of parameters, it is therefore possible to quantify the fraction of events for which 
the degeneracy can be broken, regardless of the phase of the orbit. Furthermore, 
because r does not vary with ip, for any given event the fractional errors are also 
constant throughout the orbit. 

The situation is considerably more complicated for bulge observations. For 
this case, a ~ 90°. From equation (2.2), we find that the values of r and u±_ 
vary as a function of -0- Therefore the coefficients of the constraint equation (2.5), 
and hence the ability of the constraint to distinguish between the solutions, are a 



12 



strong function of the phase of the orbit. For example, consider Baade's window, 
for which a = 84°. For ip = or 180°, we can approximate the constraint equation 
as 



whereas for i/j = 90° or 270°, 



We can see from these relations the effect of t/j on the abihty of the constraint to 
distinguish between the degenerate values of A/3. For a specified uncertainty in 
the value of Ao;, the uncertainty in A/3 is 100 times less when '0 = or 180°, 
than when ^ = 90° or 270°. Thus we can expect that, in general, it is easier 
to distinguish between the degenerate solutions when = or 180° than when 
ip = 90° or 270°. For intermediate values of ijj, we can interpolate the behavior of 
the constraint accordingly. 

For bulge observations, the intrinsic errors are also strongly dependent on ■0. 
From equation (2.2) we see that, for observations near a — 90°, the projected 
Earth-satellite separation, r, varies approximately as sinip. Thus, from equation 

(1.3), for ~ or 180°, the value of Ax is comparatively small, and therefore the 
fractional errors arc comparatively large, whereas for ip ~ 90° or 270°, the value of 
Ax is comparatively large, and thus the fractional errors are comparatively small. 
As before, we can interpolate the behavior of the fractional errors for intermediate 
values of ■0. 

The effects of 'ip on the constraint condition and fractional errors are in com- 
petition with each other. We therefore expect the DBF, which is a combination 
of both of these considerations, to have four minima: when the fractional errors 
dominate, at '0 = or 180°, and when the constraint on A/3 is poor, at ■0 ~ 90° 
or 270°. We expect the maximum values of the DBF to occur where the fractional 
errors are not dominant, yet the constraint still distinguishes between the degen- 
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erate solutions reasonably well. In the next section, we review these qualitative 
arguments in a more precise form. 

3. Error Analysis 

In this section we give an overview of the mathematical analysis underlying the 
Monte Carlo code. For a more complete description, we refer the reader to § 2 of 
Boutreux & Gould (1996); much of the analysis given therein is entirely applicable 
to the present considerations. The reader who is primarily interested in the results 
can skip this section entirely. 

The flux F{t) from a microlensed star is a function of flve parameters, = 
{to, P, uj, Fq, B), with the form, 

F{t; to, (3, u, Fq, B) = FoA[x{t; to, (3, u)\ + B, (3.1) 

where Fq is the unlensed flux, and B is the light from any additional unlensed 
sources that are not resolved. The flux from the satellite can also be described 
in this manner, by flve additional parameters a[ — {tQ, P' , cu' , Fq, B') . Using mea- 
surements simulated by our Monte Carlo method, we then flt the parameters by 
minimizing x^- By differentiating in the neighborhood of the solution, one can 
obtain Cy , the covariance matrix of the for both the Earth and the satellite. 

We assume that the Earth and sateUite have the same fllters, implying that 
Fo = Fq and B — B'. Each of these constraints can be written ctjaj = 0, 

where are now the ten parameters for the Earth and the satellite, and ai are 
the constraint coefficients. The constraint coefficients for Fq = Fq are ccj = 
(0, 0, 0, 1, 0, 0, 0, 0, —1, 0). We impose these constraints on the 10 x 10 covariance 
matrix Cy, forming a new matrix Cy, 

Cij — Cij 1Q . (^o.zj 

Neglecting the relative Earth-sateUite motion, we find that there is a two-fold 
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degeneracy in the magnitude \Af3\ = Af3±, and another two-fold degeneracy in the 
sign. Thus there are four possible solutions, ±Aj3±. In our simulations, we know 
the true solution, which for the illustrations given here we take to be +A/3_ or 
( — h) for short. However, we initially assume that only the magnitudes of (3 and 
P' are known, thus making all four possible values of A/3 equally probable. We 
use the constrained covariance matrix, Cy, to form a new covariance matrix of the 
quantities Oj = {At, A/3, Acu); we call this new matrix c""*". Due to the fact that 



there are four possible values of A/3, we form four different matrices c^-^. We then 
test to see if the solutions are distinguishable from the true solution, using the 
constraint imposed by the relative motion of the Earth and the satellite, equations 
(2.3) - (2.5). 

The general form of can be written as 

xu - ■ (3.3) 

Z^m,n=l (^mnOimOin 

where the ckj's are the coefficients from the constraint condition, and 5/3 is the 
difference between the true and trial values of A/3. Because is a strong function 
of the phase of the orbit, it is not possible to give iii a simple form. It is 
illuminating, however, to examine the form of the values of ^ for which it is 

at maximum and minimum. For ip = (90°, 270°), the constraint coefficients are: 

ctj = (0, cos a, 1) (3.4) 



and for t/j = (0,180°), 



ai — {0,fl^, cos a). (3.5) 



We define q = (ai^/apte^^), where = Cgg'^tg, and ap = C22 ■ Using these ex- 
pressions, and assuming (as is generally the case), that the off-diagonal covariances 
(cy, % 7^ j) are small, the equation for can be written for each case as. 
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,2 ,,2/nno o^noN 



Xmin = X'(90°, 270°) 



l + (g/cosa)2 ^^^^ 
y2 = v2(0 180°) = ^^/^^ 

For typical events, ~ 6 for the ±A/3_ solutions; this value is roughly independent 
of ■0- 

Note that Xmin/ Xmax ~ (cosa/g)^ (for typical time scales). Thus for obser- 
vations toward the bulge, a ~ 84°, can vary by a factor of 10^ — 10^ over the 
course of the orbit. 

We calculate the for all four solutions, and consider only those solutions for 
which < 9. We then calculate the intrinsic error for these allowed solutions. 
The fractional error in the scalar solution can be approximated by {{5Ax)'^)^^'^/Ax, 
where again Ax = (a;Ai,A/?) — (01,02), and Ax — |Ax|. Using the analog of 
equation (3.2), 

+ _ + Eti,kcumcjk(yk 

we apply the constraint (2.3) to c ^ to form c ^, where ai is now given by equation 
(2.5), and Oj = (At, A/3, Acu). In terms of c^j^, the fractional error becomes 

((5Ax)2)V2_(E?,,=iC-.+ 0,0,-)V2 



Ax (oi)2 + (02)2 



(3.8) 



Prom this we can see that as the quantities (01,02) become small, the fractional 
error in the solutions rise. Recall that for ~ 0, the projected Earth-sateUite 
separation is small (relative to the true separation), and therefore Ax is small. 
Thus the intrinsic errors in the solutions become more important near — 0. 

We approximate the fractional error in the vector solution by (((5Ax)2)^/2/Ax. 
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This can be written as, 



((5Ax) 



2\l/2 



Ax 



+ C. 



-+- 



22 



(ai)2 + (a2)2 



(3.9) 



4. Bulge and Disk Models 

In order to calculate the DBF using the Monte Carlo simulation, we adopt 
appropriate models for the velocity and density distributions of the disk and bulge 
of the Milky Way. 

For the density distributions, we choose a Bahcall model for the disk. For the 
bulge, we adopt a 'revised' COBE model. This model uses a barred bulge density 
distribution for the outer parts of the bulge (r > 0.7kpc) and a Kent model for 
the inner part of the bulge. For a more detailed description, see Table 3 from Han 
& Gould (1995). 

We assume the velocity distributions to be Gaussian, with mean velocity v, and 
velocity dispersion a. For the disk, we use {vx, Vy, Vx) — (0, 220, 0); with dispersions 
of {ax,cry,az) — (40, 30, 20) km s"-*^. For the bulge, {vx,Vy,Vx) — (0,0,0), while 
((7a;, (7y, (7^) = (110, 82.5, 66.3) kms~^. The motion of the Sun is also taken into 
account, and is {vx,Vy,Vz) = (9,231, 16)kms~-'^ (Mihalas & Binney 1981). 

There are two distinct types of events that are seen towards the bulge. These 
are bulge-bulge events, where a bulge star is being lensed by bulge lens, and disk- 
bulge events, where a bulge star is being lensed by a disk lens. These two situations 
produce dissimilar events, and we therefore consider them separately. 

For disk-bulge events, the position of the source star, d^s is confined to a rel- 
atively narrow range as compared to the range of possible values for rfoi- Thus 
we can approximate the location of the source stars as being fixed at the galacto- 
centric radius, Rq = Skpc. The frequency of lensing events scales as p{r\)rcV oc 
p{r\)[z{l — z)\^/'^v, where v is the transverse velocity of the lens, r\ is the dis- 
tance of the lens from the galactic center (for fines of sight near the galactic center. 
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ri ~ dis), and p(ri) is the mass density at ri. Thus to simulate an event, we draw the 
distance to the lens randomly from the probability distribution p{rij[z{l — z)]^^'^. 
The velocities of the lens and the source are then drawn randomly according to 
the parameters described above, and the transverse velocity calculated. We weight 
the DBF by the transverse speed in order to include its effect on the event rate. 

For bulge-bulge events, the frequency of lensing events depends on the mass 
density at the position of the source, the mass density at the position of the lens, 
the size of the Einstein ring, and the transverse velocity of the lens. In addition, 
as the distance from the observer increases, the volume element along the line of 
sight also increases, and therefore so does the total number of stars. Thus the 
event rate scales as p{r s) p{ri)rQd'^^v . Here is the distance of the source from 
the galactic center, and d^g provides the contribution due to the increasing volume 
element. Thus, the source-lens pair position (dos, c^ol) is drawn randomly from 
the distribution p{rs)p{ri)[z{l — z)Y/'^d'l^. The resulting DBF is weighted by the 
transverse velocity, as above. 



5. Monte Carlo 

For the simulation, we assume that the observations begin when the source 
enters the Einstein ring radius as seen from the Earth. The observations continue 
until the source-lens separation is > 3re as seen from both the Earth and the 
satellite. The impact parameter, is either fixed, or chosen randomly in the 
range < /3 < 0.7. This range of (3 is consistent with that of observed microlensing 
events. Those events for which the time of maximum magnification for the satellite 
has already passed when the observations begin, that is, (1 — 0^)^^'^ < a; At, are 
considered unresolvable. Those events for which /3' > 1.5 are also considered 
unresolvable. 

The four possible solutions are considered separately, and the calculated 
for each. Those solutions for which > 9 are distinguishable from the true 
solution, and are therefore unallowed. We then determine if the allowed solutions 
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are sufficiently similar that it is unimportant to break the degeneracy, and if the 
fractional errors in the solutions are intolerably large. The scalar difference between 
the solutions is given by, 



S{Ax) = \AxU,^^^ 



- \Ax 



true' 



(5.1) 



whereas the vector difference is given by. 



5(Ax) = lAXaiiowed - 



-true I ■ 



(5.2) 



If the difference between all allowed solutions and the true solution is less than 20% 
of the true solution, and if the fractional error in Ax of each allowed solution is less 
than 10%, we consider that the degeneracy is broken. For each set of parameters, 
the Monte Carlo simulation was run 5000 times, resulting in errors of < 1%. 



Because the bulge is located at echptic angle a ~ 90°, there will be times of 
the year that observations will be impossible, due to the position of the Sun. To 
quantify this effect, we define an allowed range of the ecliptic longitude. A, for the 
Earth and the satellite. When the position of the Earth or the satellite is within 
this range, observations are considered impossible. The full width of this range is 
chosen to be 90° (3 months) for the Earth, and 60° (2 months) for the satellite. In 
our simulations, we consider lines of sight towards the bulge, for which A = 270°. 
The forbidden ranges for the Earth and satellite are then. 



Note that A = is the vernal equinox. The time of the year when either the Earth 
or the satellite is within the forbidden range is dependent on the Earth-sateUite 



6. Observational Considerations 



90° - 45° <Ae< 90° + 45 



90° -30° < A^ < 90° + 30°. 



(6.1) 
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separation, dsat- For dsat < 75°, the forbidden ranges of the Earth and satelhte 
overlap. Conversely, for rfgat > 75°, the forbidden ranges are separated by a gap of 
A A = dsat — 75°. Thus values of (igat > 75° are somewhat less desirable due to the 
fact that the times of year when observations are possible are separated into two 
ranges, thereby increasing the probability that the observations of an event will be 
foreshortened due to the Earth or the satellite moving behind the Sun. We will 
consider this issue more fully in § 7.1, when we discuss the optimal Earth-satellite 
separation. 

7. Results 

There are eight parameters that affect the degeneracy breaking fraction; «, the 
angle between the source and the south ecliptic pole; ip, which is related to the 
phase of the orbit; aE and as, the photometric precisions of the Earth and satellite 
observations, respectively; P, the impact parameter as seen from the Earth; dsat, 
the Earth-sateUite separation distance; M, the mass of the lens; Nmeas, the number 
of measurements per Einstein ring crossing time. 

For all of the parameters (neglecting for the moment), we choose typical val- 
ues that we consider to be realistic. These values are used in all simulations unless 
varied explicitly. A value of Nmeas = 20 corresponds to frequency of measurements 
in the range of 0.5 to 4.0 days for typical events. We choose {aE, o's) = (1%, 2%) 
as our photometric precision. The photometric precision obviously depends on the 
magnitude of the source star; these are typical values of photometric precision for 
giant sources. We choose M — O.SMq as a representative mass of lenses for ob- 
servations towards the bulge. The default Earth-sateUite separation is dgat — 60° 
(2 months). The maximum impact parameter is set at (3max — 0.7. For Baade's 
window, a = 84°. 

Observations towards the bulge are unique in that they are a strong function 
of the time of year, as discussed above. Due to this intrinsic dependence of the 
DBF on ■0, we include as a free parameter in all further calculations. 
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In order to elucidate the effect of tlie time of tlie year on tlie DBF, we present 
in Figure 2 tlie DBF against ip for our fiducial parameters. Included in Figure 2 
are the results for both the disk-bulge and bulge-bulge cases. Also included is the 
DBF against ip for disk-bulge events assuming a fractional difference threshold of 
5% (see § 2.3). Within the hmits of the Monte Carlo errors, we find that in the 
range < ip < 180° the curve is equivalent to that for the range 180° < tp < 360°. 
We utilize this symmetry by restricting the range of ip in further calculations to 
< tp < 180°. There is, however, a shght asymmetry around tp — 90° (and 
similarly around ip = 270°). For the speed degeneracy, this asymmetry makes no 
significant contribution, and we disregard it. We will consider this effect more 
fully in § 8, where we discuss the velocity degeneracy. For ip 0, 180°, we see from 
equation (2.2) that the projected Earth-satellite velocity is at a maximum, and 
as a result the constraint given in equation (2.1) is very effective in distinguishing 
between the possible solutions. The projected Earth-satellite separation is at a 
minimum, however, driving up the fractional errors in the solutions. Thus near ip ~ 
0, 180°, the DBF is dominated by the intrinsic errors. As ip increases, the projected 
Earth-sateUite separation increases, and thus the intrinsic errors decrease. The 
projected Earth-satellite velocity decreases as ip increases, and as a result the 
constraint becomes less effective in distinguishing between the degenerate solutions. 
For Ip ~ 20° the DBF is at a maximum; the constraint still distinguishes between 
the solutions for a majority of the events, while the intrinsic errors are typically 
< 10%. When ip — 90°, the projected Earth-satellite separation is equal to the 
true Earth-satellite separation, and the intrinsic errors are typically < 5%. The 
projected Earth-satellite velocity is at a minimum at = 90°, and the constraint 
distinguishes poorly among the solutions (as compared to ~ 0). Thus, a,t ip — 
90°, the DBF is dominated by the difficulty in distinguishing between the individual 
solutions. The analysis of the DBF curve in the range < ip < 90° is identical to 
that for the rest of the orbit. 

The analysis of bulge-bulge events is exactly the same as for disk-bulge events. 
We can understand how disk-bulge and bulge-bulge events differ by considering the 
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typical size of the projected Einstein ring in each case. For disk-bulge events, the 
median value for the distance to lens is rfoi ~ 5 kpc, and the distance to the source 
is dos = 8 kpc, thus, from equations (1.2) and (1.4), we find that, for disk lenses, 
fe ~ 6.OAU(M/O.3M0)V2. For bulge-bulge events, doi ~ 7kpc, dos ~ 9 kpc, 
so that fe ~ 9.OAU(M/O.3M0)V2. Therefore the typical size of the projected 
Einstein rings will be larger for bulge-bulge events by a factor of ~ 1.5. This will 
affect the value of the DBF in two distinct ways. For phases of the orbit such that 
■0 ~ 0, where intrinsic errors dominate, the DBF for bulge-bulge events will be 
smaller than that for disk-bulge events, due to the larger typical values of fe (and 
therefore larger fractional errors). For phases of the orbit such that ip ~ 90°, we 
expect that the DBF for bulge-bulge events will be, in general, higher than the 
corresponding DBF for disk-bulge events, due to the fact that the degenerate Aj3 
solutions are better separated in the case of the bulge-bulge events than those for 
disk-bulge events. 

For disk-bulge events, we find that for our fiducial parameters, it is possible 
to resolve > 87% of the events during most times of the orbit, and > 70% of the 
events over the entire orbit. For bulge- bulge events, we find values of > 90% for 
most of the orbit and > 60% for the entire orbit. 

Figure 3 represents the DBF as a function of i/j and /3 for disk-bulge events. 
For this figure, [3 is not randomly chosen, but rather held fixed for each point in 
the curve. Contrary to the results for the LMC (Boutreux & Gould 1996), where 
it was found that events with /? > 0.7 are much harder to resolve, there is no one 
value of (3 for which the DBF drops below a certain value. In fact, for — 90°, 
even those events for which P ~ 0.9 are resolvable 80% of the time, whereas less 
than 20% of these events are resolvable when t/j — 10°. In other simulations, we 
choose P randomly between < P < 0.7 because the majority of events seen so far 
have values of P that are in this range, but we emphasize that, under the proper 
conditions, even those events that have /5 > 0.7 may still be resolvable. 

Figure 4 represents the DBF as a function of t/j and M for disk-bulge events. 
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Figure 2. Degeneracy Breaking Fraction (DBF) as a function of ijj for disk-bulge events (solid bold curve) 
and bulge-bulge events (dashed curve). All other parameters are fixed at Uobs = 20, [oe-^'^s) = (1%,2%), 
M = O.3M0, dsat = 60°, (5max = 0.7, and a = 84°. The DBF is the fraction of events for which the allowed 
solutions differ by < 20% and each of these has intrinsic error < 10%. The lower solid curve is the DBF as 
a function of ij) when allowed solutions differ by < 5%, and each has intrinsic error < 10%. 

For the mass, we focus on the range O.OIMq — IM©; this corresponds to the best 
estimate of the range of masses for events already observed. For our fiducial value, 
M = 0.3Mq, the DBF is ~ 90% for the majority of the orbit. For disk-bulge 
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disk— disk 




Figure 3. DBF as a function of ip and /3 for disk-bulge events. The labeled bold contour represents a 
DBF of 90%. Other bold contours represent DBFs of 50% and 10%. Intermediate contours are spaced at 
intervals of 10%. All other parameters are held fixed at our fiducial parameters, as in Fig. 2. 

events, the projected Einstein ring radius is fe ~ 6.OAU(M/O.3M0)V2. Consider 
first phases of the orbit for which ip ^ 0, where the fractional errors dominate 
the behavior of the DBF. For these values of i/j (and R ~ lAU), the projected 
Earth-sateUite separation is r ~ O.lAU. The fractional errors become important 
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when Ax <^ 1. A quick calculation shows that, for masses greater than M > 
O.IMq, the magnitude of the vector displacement is Ax = r/fg < 0.1, and thus 
the fractional errors begin to dominate, thereby lowering the DBF. For masses 
M < 0.03Mq, Ax = r/fe ~ 1, and thus the fractional errors are unimportant, and 
the DBF is > 90%. Consider now phases of the orbit for which ~ 90°. For the 
range of masses considered here, Ax > 1, and thus the fractional errors are never 
large. Considering increasingly smaller masses, we find that more of the events 
are unresolvable, due to the fact that the projected Earth-sateUite separation, r, 
becomes increasingly larger than the Einstein ring radius, fe, and the event as seen 
from the satellite will more frequently fall outside of the Einstein ring. Thus the 
DBF drops rapidly from 90% to 40% from M = O.SM© to M = O.OSM©. 

Figure 5 represents the DBF as a function of -0 and B, the echptic latitude (B = 
a — 90°), for disk-bulge events. For this figure, we hold the echptic longitude fixed 
at A = 270°, and vary B in the range —10° < B < 10°. In galactic coordinates, 
this corresponds to the range (Z, b) — (—2°, —5°) — (16°, 4°). For Baade's window, 
B = -6°, the DBF is greater than 90% for the entire orbit. Note that for B = 0, 
it is impossible to break the degeneracy. This is due to the fact that, at B = 0, the 
constraint given in equation (2.1) becomes X^i^i '^i'^i — + ^©t^At cot = 0, 
and thus does not constrain the value of Aj3 at all. Similarly, for small values of 
|B|, the value of Af3 is poorly constrained. Thus, there is range of ecliptic latitude, 
|B| < 1°, where the DBF is less than 40% for a majority of the orbit; this range 
should be avoided when making observations. Lines of sight where |B| > 2° are 
optimal. 

Figure 6 represents the DBF as a function of i/j and (igat, the Earth-satellite 
separation, for disk-bulge events. Also shown is the range of ip for which obser- 
vations are impossible due to the position of the Sun. This range is a function 
of (igat; thus by considering simultaneously how this range and the DBF vary as 
function of dgat, we can find an optimal Earth-satellite separation. For dgax < 30°, 
the range where observations are impossible is smallest (~ 90°), but the DBF is 
never greater than 80%. For dgat — 60°, the observationally forbidden range is 
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Figure 4. DBF as a function of and log M/M© for disk-bulge events. Contours are the same as those 
in Fig. 3. 

only slightly larger (~ 130°), and the DBF is greater than 90% for the majority 
of the orbit. For dsax > 75°, the forbidden ranges for the Earth and the satellite 
begin to diverge, and the times of the orbit where observations are possible become 
separated. Thus, although the DBF continues to improve as dgat increases, there is 
an increased likehhood for these ranges that an observation will be foreshortened 
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Figure 5. DBF as a function of ip and B for disk-bulge events. Contours are the same as those in Fig. 

due to the Earth or the sateUite moving behind the Sun. We consider 20° — 75° to 
be the optimal range for the Earth-sateUite separation. For dgat < 20°, the DBF 
is unacceptably low. 

We also considered the dependence of the DBF on the parameters M, B, 
and dsat for bulge-bulge events. We find the results to be almost identical to those 
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Figure 6. DBF as a function of ij) and rfgat for disk- bulge events. The y-shaped area of the graph between 
the heavy hnes represents those range of for each value of dsat where observations are impossible due to 
the position of the Sun. The other contours are the same as those in Fig. 3. 



for disk-bulge events, with variations of at most ~ 10%, when the phase of the 
orbit is ~ or 180°, but typically the variations are much smaller. For reasons 
discussed above, the DBF is, in general, somewhat higher when -0 ~ 90°, and 
somewhat lower when ~ for bulge-bulge events. As with disk-bulge events, we 
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find that there is no natural cutoff value for for our fiducial value of M = O.SMq, 
the DBF is ~ 90% for the majority of the orbit; for Baade's window, B = —6°, 
the DBF is greater than 90% for the majority of the orbit; and the optimal range 
of o?sat is 20° — 75°. Overall, the differences between bulge-bulge and disk- bulge 
events are minor, and will not significantly affect the observations. 

8. The Velocity Degeneracy 

In the previous section, we discussed the observational requirements of a satel- 
lite to break the two-fold degeneracy in the projected speed, v. While measuring 
the projected speed is helpful in distinguishing between component populations, 
a measurement of the projected velocity, v, would provide a more conclusive dis- 
tinction. 

In this section we present the results for breaking the four-fold velocity degen- 
eracy, adopting the fiducial parameters used in the previous section: Nmeas — 20, 
((TE,(Ts) = (1%,2%), M = 0.3Mq, dsat = 60° (2 months), (3max = 0.7, and 
a — 84°. We calculate the difference between the allowed solutions using equation 
(5.2), and the fractional error in each allowed solution using equation (3.9). As 
in the previous section, if the intrinsic errors in all allowed solutions is < 10%, 
and if the difference between the allowed solutions is < 20% of the true solution, 
we consider the degeneracy to be broken. Here we present the results for both 
disk-bulge and bulge-bulge type events. 

Figure 7 represents the vector DBF as a function of t/j for disk-bulge events 
and bulge-bulge events. The scalar DBF is also shown for comparison. We find 
that for disk-bulge events the DBF is greater than 40% for the entire orbit, 
and greater than ~ 55% for about half the orbit. We find that for bulge-bulge 
events the vector DBF is ^ 30% over the entire orbit, and ^ 40% over half the 
orbit. In general, the analysis of the vector DBF curve is the same as that for 
scalar DBF curve: there are minima at = and 180° where the intrinsic errors 
dominate the DBF, and a minimum at = 90° where the ability of the constraint 
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Figure 7. Vector DBF as a function of ip for disk-bulge events (bold solid curve) and for bulge-bulge 
events (bold dashed curve). Also shown is the corresponding scalar DBF. All other parameters are fixed 
at nobs = 20, {(7E,as) = (1%,2%), M = O.SMq, 4at = 60°, Pmax = 0.7, and a = 84°. 



to distinguish between the degenerate solutions dominates the DBF. There are 
two unique features of the vector DBF, however. The first is that the minimum 
value of the DBF occurs at ~ 90°, where the limiting factor is the ability to 
distinguish between the solutions. In the scalar case, distinguishing between the 
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±A/?_ solution is unimportant, due to the fact that these solutions have essentially 
the same projected speed. To break the vector degeneracy, however, one must 
distinguish between the ±A/9_ solutions, which is often difficult when ip ~ 90°, 
where the constraint is the least effective. We also find a significant asymmetry in 
the vector DBF curve around ^0 = 90°, the largest asymmetry being a difference 
of ~ 15% between the DBF at = 15° and 165°. This asymmetry is caused by 
fact that lenses in the disk have a preferred direction of motion, and, to a lesser 
extent, from the motion of the Sun around the galactic center. For a given event, 
A/3 oc vx, where v±_ is the component of the transverse velocity perpendicular to the 
projected Earth-satellite separation, r. As changes, the direction perpendicular 
to r changes with respect to the preferred direction of disk. Thus a preferred 
direction for v will translate into an asymmetry in the distribution of values of A/3 
around 'ip = 90°. Specifically, the values of A/3 are typically smaller when ip < 90° 
than when -0 > 90°. Therefore, because the ±A/3_ solutions are, on average, more 
closely spaced when t/j < 90°, the constraint is less likely to distinguish between 
them, and the DBF is smaller. This asymmetry does not arise in the scalar DBF 
curve because it is not necessary to distinguish between the ±A/3_ solutions. Note 
the smaller asymmetry around = 90° for bugle-bulge events. This is due to the 
fact that, in bulge-bulge events, there is no preferred direction of motion of the 
lens as there is in disk-bulge events. 
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APPENDIX Generalized Satellite Orbits 



For our Monte Carlo simulations, we assumed that the satellite is in an Earth- 
like orbit, and that both the orbit of the Earth and the orbit of the satellite are 
circular. In reality, a satellite launched from the Earth will not be in an Earth- 
like orbit, but will actually be in an elliptical orbit. This implies that the true 
velocity of the satellite will differ from the assumed velocity for each value of the 
Earth-satellite separation. In this appendix, we show that, for a reasonably small 
deviation from a circular orbit, the difference between the true and assumed Earth- 
sateUite velocity will be small, and that therefore our results are sufficiently generic 
to be apphcable to a large class of orbits. 

When launched from the Earth, the satellite will have some velocity with re- 
spect to the Earth. We assume the direction of the satellite's velocity is parallel to 
that of the Earth (i.e. the orbit of the satellite is in the plane of the Earth's orbit, 
and the satellite's perihelion is at the distance of the Earth). That is, Vgat — v^+5v. 
We assume that the orbit of the Earth is circular, therefore a non-zero value oi Sv 
implies that the satellite will have an elliptical orbit with a period different from 
that of the Earth's, and it will drift away from the Earth over time. Using Kepler's 
laws, one finds that the difference between the period of the Earth and the satellite 
is given by, 

'-^ = (A.1) 

Pa 

In § 8.2, we determined that the optimal Earth-satellite separation was in the range 
20° — 75°, or I month — 2 ^ months. We would like the Earth-satellite separation 
to vary in this range over the course of ~ 3 years. For definiteness, we choose the 
Earth-satellite separation to be 3 months after 3 years, i.e. 5P/Pq ~ 1/12. Thus, 
from equation (A.l), and using = 30kms~-^ we find 5v ~ O.Skms"-^. The orbit 
is nearly circular, with an eccentricity of e = 0.06. 

Using the parameters of the orbit, we can now calculate the difference between 
the true value of the Earth-satellite velocity and the value used in our simulation 
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Figure 8. Earth-satcUitc velocity as a function of time from launch. Shown are the Earth-satellite velocity 
assuming a circular orbit (dotted line), and the real Earth-satellite velocity (solid line). 



(which assumed a circular, Earth-hke orbit for the satelhte). We set t = as the 
time of launch, and at each point in time, we calculate the true Earth-satellite 
separation and velocity. We then calculate what the value of the Earth-satellite 
velocity would be if the satellite were in an Earth-like orbit, separated from the 
Earth by an amount equal to the true Earth-sateUite separation. This corresponds 
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to the velocity used in our Monte Carlo code. Figure 8 shows the Earth-satellite 
velocity for both cases in the range < t < 3 years. Focusing on t > 0.5 year, we 
find that the true velocity is systematically smaller than the assumed velocity by 
~ 6%, but varies from being ~ 18% smaller to ~ 6% larger. 

Here we have calculated the difference between the true Earth-satellite velocity 
and that assuming an Earth-like orbit for the satellite. For the purpose of degener- 
acy breaking, however, it is the projected Earth-satellite velocity that is important. 
In the case of an elliptical orbit, the magnitude of the projected Earth-satellite ve- 
locity is dependent on the time of year in which the satellite is launched. Since we 
do not know a priori when this will be, we cannot find precisely the difference be- 
tween the true and assumed projected Earth-satellite velocity. This is not strictly 
necessary, however, since the difference between the projected velocities will be of 
the same order as the difference between the space velocities. Since the ability 
to break the degeneracy is roughly proportional to the projected Earth-sateUite 
velocity, we conclude that the error introduced in our simulation by assuming an 
Earth-like orbit for the satellite is of order the error in the Earth-satellite velocity, 
~ 6%. 

We conclude that the results obtained in our simulations are applicable to 
likely sateUite orbits. If, however, one would hke to consider satellite orbits that 
are either highly inclined to the Earth's orbit, or highly eccentric, then our results 
are not applicable. These situations will need to be analyzed individually. 
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